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ABSTRACT 


Topaz II is a Russian made nuclear powered satellite purchased by the United States. This 
thesis provides a vibration analysis of the power unit of Topaz II. Several experimental vibration 
tests were conducted on the actual power unit by Sandia National Laboratories, in Albuquerque, 
New Mexico. Utilizing Structural Dynamics Research Corporation’s I-DEAS Masters Series 
Software, a finite element model was made of the power unit. Using the I-DEAS software, this 
model was tested using the same vibration test values as those performed on the actual power 
unit. The lateral and axial, experimental and theoretical, sine tests were then analyzed. The 
purpose of this analysis was to determine how accurate a representation the model was of the 
actual Topaz unit. This was accomplished through comparison of natural frequencies obtained 
experimentally with those predicted from I-DEAS testing of the finite element model. The 
frequency range of interest was from 0 Hz to 200 Hz. When solved for system modes, the model 
predicted all natural frequencies within 4.6 Hz of those found experimentally. The average 
variation between theoretically predicted and experimentally determined natural frequencies was 


1.8 Hz. 
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I. INTRODUCTION 
A. ACQUISITION OF TOPAZ I 

The United States conducts many of its military missions in space through the use of 
satellites. The Department of Defense has a long term requirement to develop follow-on power 
technologies to support these missions. Having the capability to produce a large quantity of 
electrical power on orbit would allow a greater spectrum of missions to be performed, such as 
environmental monitoring and using high power active sensors for surveillance. To determine the 
most effective means of supplying this power, a series of Air Force and Department of Defense 
mission and integration studies were begun in the late 1980's. The conclusion drawn was that 
compact, environmentally safe nuclear power was the best option. In the spring of 1992, a 
Department of Defense study considered the option of importing the Russian made TOPAZ II 
space nuclear reactor technology. It was decided that this technology would greatly enhance 
development of U.S. space nuclear reactors. [Ref. 1] 

Procurement of TOPAZ II was accomplished by the Strategic Defense Initiative 
Organization (SDIO), now known as the Ballistic Missile Defense Organization (BMDO). They 
have also been responsible for the transfer of space nuclear reactor technology between the 
United States and Russia. This started in December of 1991, when SDIO began the investigation 
of the possibility of a U.S. launched Russian TOPAZ II nuclear power system. At this time, the 
intended destination of TOPAZ II was in the Nuclear Electric Propulsion (NEP) Space Test 


Mission. [Ref. |] 
The foremost consideration in launching a TOPAZ I satellite was safety. The 


requirement was for a safe launch, within constraints of the given budget. To find out if this were 


possible, SDIO subpoenaed a group of twelve engineers and scientists from Sandia National 
Laboratories (SNL), the Air Force Phillips Laboratory (PL), Los Alamos National Laboratory 
(LANL), the University of New Mexico (UNM), and Advanced Sciences, Inc. (ASI). The 
conclusion of this comprehensive preliminary study, completed in 1992, was that the TOPAZ II 
could be launched safely and within budget constraints. 

In the spring of 1992, two TOPAZ II reactors were delivered to the New Mexico 
Engineering Research Institute (NMERI). The fuel elements were not included in the package. 
Also obtained by NMERI was a complete satellite qualification testing facility. Under the 
direction of the Phillips Laboratory, non-nuclear testing of TOPAZ II was begun as part of the 
TOPAZ Technology Program. 

The TOPAZ Technology Program, led by the BMDO, is a three phase cooperative effort 
with Russia. The three phases are comprised of joint evaluation, joint testing, and joint 
development. Each phase allows an opportunity for both the U.S. and Russia to evaluate the 
program and determine future direction. The United States stands to gain possession of 


thermionic reactor technology for use by U.S. industry. [Ref. 1] 


B. MISSION DESCRIPTION 

Mission goals of the NEP Space Test Program include use of the NEP Space Test 
Spacecraft to test an assortment of electric propulsion thrusters and conduct additional scientific 
experiments that would be compatible with its other mission requirements. Two specific elements 
of Russian space technology will also be tested and evaluated, those being the TOPAZ II nuclear 


reactor and the SPT-100 electric propulsion thruster. One of the programs primary goals is to 


evaluate the NEP in orbital transfer applications and measurement of the spacecraft environment 


produced by NEP system operation. 


The reactor system of the NEP spacecraft, shown in Figure 1, is separated from the main 


section of the spacecraft by an extendable boom. This arrangement is necessary to maintain the 


radiation dosage received by other spacecraft components to acceptable levels. [Ref. 1} 


C. TOPAZ Il POWER SYSTEM DESCRIPTION 


The power system of TOPAZ II consists of a reactor power system, which generates 


electricity from nuclear heat through the use of in-core thermionic conversion units. This system 


was designed by Russian engineers to meet the following specifications: [Ref. 1] 


L. 


Under no circumstances should the reactor system operate before reaching orbit 
altitude, 

The system shelf life must be 10 or more years 

The coolant shall not freeze before operation 

The system must be capable of providing 6 kWe at 27 volts for a lifetime of 3 years 


with a reliability of 0.95. 


. The mass of the power system must not exceed 1061 kilograms, excluding the 


controller and supply system masses. 


’ The main subsystem of the TOPAZ II power system are the reactor subsystem, the radiation 


shield, the primary coolant loop, the cesium supply system, the instrumentation and control 


system, gas systems, thermal cover, and the primary power system structure. (see Figure 2) 
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Figure t. TOPAZ [Ref. 1]. 


The nuclear reactor is a small zirconium hydride-moderated, epi-thermal reactor with in- 
core thermionic converters. The reactor incorporates 37 in-core single cell thermionic fuel 
elements (TFEs). Prior to fuel being loaded into the TFEs, electric heaters can be loaded into the 
TFE cavity to simulate heat generated from the reactor. The advantage of using heaters is the 
capability to perform non-nuclear testing of the thermionic converters and the complete power 
system. With heaters, this can be accomplished at conditions near what the nominal operating 
conditions would be. 

The purpose of the radiation shield is to attenuate neutron and gamma radiation. Located 
behind the reactor, it is a truncated cone composed of a stainless steel shell containing lithium 
hydride (LiH) which is the neutron shield. The radiation shield is attached by support legs to the 
lower end of the reactor. 

The reactor coolant loop is comprised of a single EM pump, stainless steel piping, sodium- 
potassium (NaK) coolant, and a heat rejection radiator. From a lower plenum, the NaK passes 
through the reactor core where it is heated, to an upper plenum to two stainless steel pipes to the 
radiator, the NaK flows through two coolant pipes dividing into three pipes each, prior to 
reaching the EM pump. The pump then sends the Nak back to the reactor lower plenum. 

Cesium is supplied to the TFE interelectrode gap via the cesium supply system. This is 
necessary to increase the efficiency of the thermionic converters through suppression of the space 


charge that occurs near its emitters. 
The Instrumentation and Control system is required for the monitoring, controlling, and 
telemetering of the power system status. Its primary duties are the following: 1) to start the 


power system, 2) maintenance of the system operation under nominal operating conditions, 3) 
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Figure 2. Primary Power System [Ref. 1]. 


stabilization of the voltage received by the payload 4) performance of the commands from the 
ground control station, 5) shutdown of the power system, 6) maintain safety control during land- 
based operations, 7) telemetering performance data to the ground, 8) shunting of excess 
electrical power to ballast resistors, and 9) to keep the storage battery charged. [Ref. 1] 

There are four gas systems on board TOPAZ II. These are the CO,/He moderator cover 
gas, the He TFE-gap gas, the Ar/He gas in the volume accumulator and the helium gas in the 
radiation shield. The gas pressure is monitored during filling, electric ground testing, and 
throughout the first few thousand hours of space operation by pressure gauges connected to each 
of these gas systems. 

The purpose of the thermal cover is to reduce the heat loss from the power system on the 
launch pad and during the flight to orbit. If the TOPAZ II coolant were allowed to freeze, 
coolant circulation would cease, and the system would no longer operate properly. 

The primary power system structure consists of the reactor, the shield, and the frame. 
These are the three main structural members of the power system. The frame supports the 
radiator, the gas bottles, coolant piping, and the volume accumulator. It is attached to the shield, 
and is also used to attach the spacecraft boom. [Ref. 1] 

A more indepth analysis of the nuclear power system was done by LT Steve Benke and 
LT Rich Venable. Their theses may be obtained through the Naval Postgraduate School, 


Monterey, CA. [Ref. 2]. 


Il. TOPAZ Il - THE MODEL 
A. DEVELOPMENT / PURPOSE OF TOPAZ II FINITE ELEMENT MODEL 
The Finte Element Model of TOPAZ II was developed by J.M. Lacy, of Idaho National 
Laboratories [Ref. 3]. One of the first problems encountered was that of which system of units 
to use. American engineers tend to favor the SAE system whereas the Russians use their version 
of the SI system. Mr. Lacy elected to use the Russian modified metric system. The following 


units conversion table may prove useful. [Ref. 3] 


Unit Dimension Equivalencies 

length 1mm 0.03937 in 0.00329 ft 
mass lkg 2.2 Ibm 0.00571 Ib-s?/in 
force 1N 0.2248 Ib 0.1019 kgf 
stress 1 MPa 1 N/mm? 145.04 psi 


The finite element model was made for the reactor unit and its support frame. The frame 
has three legs whose feet, when connected by an imaginary circle, form a circle of 2440 mm in 
diameter. The upper end of the frame, to which the reactor unit is attached, has a diameter of 
680 mm. This platform is 2116 mm above the feet. (see Figure 3) On the launch pad, the power 
system and frame will weigh approximately 1043 kg. The center of gravity will be 2304 mm 
above the feet. 

Some of the TOPAZ frames are constructed of 12X18HI10T stainless steel while others 


are reported as 08X18HI10T material. These designations are interpreted below and their 


properties are listed in Table 1. [Ref. 3] 


























12 less than 0.12% Carbon 


X18 18% Chromium 
H10 10% Nickel 
T Titanium is 5 times the Carbon content, 
but less than 0.8% 
density, p 7.95 g/em%3 
modulus of elasticity, E 195.12 GPa (28.3 x 10° psi)* 


* an assumption, to correspond with ASME 
stainless steel properties 


Ultimate Yield Reduction 
Designation | Temperature Strength Strength Elongation in Area 








Table 1. TOPAZ Il Frame Properties |Ref. 3]. 


One purpose of creating the FEM, and the one addressed by this thesis, is the 
determination of the natural frequency of the reactor system. This is a major concern during the 
launch phase, due to the vibration forces that will be imparted to TOPAZ. The natural frequency 
determines how the system is going to respond to these inputs and whether or not it will 
withstand the launch loads. If the natural frequencies are accurately calculated, then the stiffness, 
mass, and mass distribution of the system are known to have been modeled correctly. Having 
accomplished this, the model can then be used reliably to predict what stresses and displacements 
will occur under actual loading conditions. 

Through testing of an actual TOPAZ unit in 1992, the Russian scientists found the first 


natural frequency of the reactor unit to be 4.5 - 5.0 Hz. This presented a problem because the 
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natural frequency of the entire system would be far lower than that of the rector unit alone. Due 
to interference from launch loads, a payload whose natural frequency is less than 4 Hz cannot be 
launched in the United States. 

When E.D. Schaefer of APL [Ref. 8] created a FEM to confirm the frequency range of 4.5 
to 5 Hz, he instead found the first natural frequency of the system to be closer to 10.5 Hz. Early 
in 1993, Lacy [Ref. 3] was tasked to study this problem independently. His results were the same 
as Schaefer's, At this point, the discrepancy explanation was thought to lie in differences between 
the reactor system and its documentation. Through hands-on measurements, this was disproved. 
Properties that could not be directly measured were confirmed with the Russian scientists. 

During the summer of 1993, Sandia National Laboratories performed a modal survey of 
an actual reactor unit. This system massed 899 kg (as opposed to 1043 kg used the previous 
FEM) and the first natural frequency was found to be 20.3 Hz. In their testing, Sandia used both 
the reactor unit and the seismic mass. When this mass was included in Lacys' FEM and the 
reactor mass reduced to 899 kg, he predicted the first natural frequency to be 21.6 Hz. 

Therefore, Sandia's modal survey supported both Lacys' and Schaefer's theoretical results of 
10 Hz. 

Lacy produced another FEM of TOPAZ in the fall of 1993 [Ref. 3]. This model weighed 
899 kg and had small refinements made to the mass distribution of the system, and to the frame 
geometry. His results continued to agree with Schaefer's results and those of Sandia National 
Laboratories. The first natural frequency of the seismic mass was only 2% below that reported by 


the modal survey, and the first lateral fixed base frequency of the system was found to be 10.8 Hz. 


(Ref. 3]. 
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The results presented in this work were derived through vibrational analysis of this 


updated FEM model. 
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I. EXPERIMENTAL SHOCK AND VIBRATION TESTING OF TOPAZ II 
A. SETUP 

Early in September, 1994, shock and vibration testing was performed on TOPAZ II at 
Sandia National Laboratories in Albuquerque, N. M. Specifically, the vibration tests performed 
were lateral and longitudinal sine vibration tests, random vibration tests, and shock tests. The sine 
tests are examined in this thesis. Analysis of the random vibration tests has been performed by 
LT Sheryl Campbell from the Naval Postgraduate School, in a companion thesis completed in 
June 1995. | 

There were twelve accelerometers mounted on the power system. The accelerometers 
were placed where the anticipated responses were deemed most critical. Their locations are given 
in Tables 2 and 3 for the lateral and longitudinal tests respectively. Also shown in this table are 
the instrument axes from which responses were read by the testing equipment. Data were not 
taken from all axes from each accelerometer during each test because the testing equipment could 
only read a maximum of 24 channels concurrently. 

There were three different coordinate systems used between the experimental set-up and 
theoretical analysis. The actual TOPAZ unit, the FEM, and the accelerometers each had their 
own coordinate directions. The relationship between the instrument and accelerometer axes is 
given in Table 4. Figure 4 shows the relationship between the FEM and TOPAZ coordinate axes. 
Table 4 does not include accelerometer numbers 11 and 12. These accelerometers were mounted 
on the cesium unit and the start-up unit of TOPAZ. In the FEM, these units were not individually 


modeled so that there was no theoretical result obtained from these Jocations. In future research, 
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Instrument Location 


Mount on bottom collector of radiator on +Y axis 


Mount on -Z leg of frame at top of radiator 





Mount on leg between +Y and +Z of frame at top of 
radiator 








Mount on reactor leg bracket most closely aligned 
with -Z axis 


Mount on reactor leg bracket most closely aligned 
with +Y axis 


Mount on reactor top plenum on -Z axis x 














Table 2. Lateral Sine Sweep Accelerometers [Ref. 6]. 


these units should be modeled and their responses evaluated. 

In this thesis, the terms lateral and longitudinal refer to the lateral (Z) and axial (X) 
directions of the actual TOPAZ unit. The position of the accelerometers on TOPAZ are depicted 
in Figures 4 and 5. Also shown are the coordinate systems of both the actual TOPAZ unit and the 
FEM, for future reference in later chapters. 

The values of the sine vibration applied to TOPAZ were chosen as common launch 
frequencies that might be imparted to the satellite. The specific values and test parameters are 


listed below. (A graphical repesentation of the sine input is given in Figure 6) [Ref. 6] 
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Instrument Location 


Mount on -Z leg of frame at base x 
Mount on leg between + Y and +Z of frame at base ae eel 


Mount on bottom collector of radiator on -Z axis 
Mount on bottom collector of radiator on +Y axis 


Mount on -Z leg of frame at top of radiator 


Mount on leg between + Y and +Z of frame at top of % x 
radiator 


Mount on reactor leg bracket most closely aligned 
with -Z axis 
Mount on reactor leg bracket most closely aligned 


with +Y axis 


Mount on reactor top plenum on -Z axis 


Mount on reactor top plenum on +Y axis 
Mount to "hard point" on cesium unit 


Mount to "start-up unit" frame 








x | 


x 














Table 3. Longitudinal Sine Sweep Accelerometer Locations [Ref. 6]. 
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SINE-VIBRATION INPUT 





Frequency(Hz) Level: g (1 min. @ axis) 
5 25 
5-8 linear increase 
8-40 1.0 
40 - 100 9 
100 - 200 8 


1. The input was ramped 

2. The test duration (sweep rate) was 0.25 octaves/minute which corresponds to a test 
duration of about 7 minutes. 

3. A standard test tolerances of +/-10% was placed on the test parameters. 

4. Three additional control accelerometers were used. Their purpose was to limit the 
response of the TOPAZ system to prevent overstressing it. These were accel- 
erometers number 2, 5, and 9. See Figures 4 and 5 for their locations. [Ref. 7 


These controls were limited to responses of 6, 9, and 12 dB respectively. [Ref. 9]. 


The sine vibration was applied to the TOPAZ unit twice, once in the longitudinal X 
direction and once in the lateral Z direction. Due to the geometry of TOPAZ, the Z direction was 


determined by the Russians to be the most critical of the two lateral axes directions for testing 


purposes. 
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On Node | Instrument Accelerometer ID 


101/C13:-Z LEG@BASE Axis: -X 
101/C14:-Z LEG@BASE Axis: -Z 
101/C15:-Z LEG@BASE Axis: -Y 

102/C16:LEG/+Y AND +Z/BASE Axis: 























105/C33:JOINT ON LEG Axis: -Z 
106/C19:LEG +Y AND +Z/RAD.TOP Axis: -X 


106/C20:LEG/+Y AND +Z/RAD.TOP Axis: 


+Z(30) 





107/C05:REAC.LEG BRACKET/-Z Axis: 
degree -X(20) 


7 79 Z -Z(2 107/C08:REAC.LEG BRAK./-Z Axis: -Z(15) 


8 75 x -X(20 | 108/CO1:REAC.LEG BRACKET/+Y Axis: 
degree -X(20) 
























108/C04:REACT.LEG BRAKET/+Y Axis: +Z 


Accel | On Node | Instrument | Topaz Accelerometer ID 


# # Axis Axis 


8 75 Z +¥ (15 








109/C09?:REACT.TOP PLENUM/-Z Axis: -X 
-109/C11:REACT.TOP PLENUM/-Z Axis: +Y 
109/C12:REAC.TOP PLENUM/-Z Axis: -Z 




















Table 4. Correlation of Coordinate Systems. 
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Figure 4. Location of Accelerometers, +XY Topaz Plane [Ref. 7]. 
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Figure 5. Location of Accelerometers, +XZ Topaz Plane [Ref. 7]. 
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Figure 6. Sine-Vibration Input. 
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B. DATA MANAGEMENT 

The data received from these tests were converted over to DOS readable ASCII format. 
Results were then manipulated by a simple MATLAB program and plotted as acceleration versus 
frequency. One such plot was obtained for each accelerometer axis from which data were taken, 
for both the lateral and longitudinal sine tests. These graphs are contained in Appendix D 

Each of these graphs has been analyzed to determine at which frequencies peaks occurred. 
These peaks represent natural frequencies of the TOPAZ system which were excited by the 
vibration input. Below 100 Hz, most peaks are easily distinguishable. Above 100 Hz, there are 
many smaller peaks appearing in the data. This was due to a combination of noise in the 
equipment and interference from excitation of lower system modes. Only prominent peaks are 
tabulated. 

One example of these graphs is given in Figure 7. Peaks are indicated by 'x' marks. The 
corresponding data aregiven in Table 5. The values were obtained by first analyzing the plot and 
then going into the converted data files to locate the actual frequency and amplitude values. One 
such plot and table were made for each accelerometer axis from which data was taken, for each 
test. The remaining graphs and tables may be found in Appendices D and E, respectively. 

The system was tested up to 200 Hz. Frequencies above this value are insignificant during 
a launch scenario and are typically determined during acoustic testing. They are also more 


difficult to accurately distinguish. 
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Amplitude, g 


Transducer: 111/C30:HARD POINT/CESIUM UNIT Axis: +Z 











fe) 2 : 3a 
10 10 10 " 
Frequency ~ Hz 


Figure 7. Example Experimental Data Plot- 


bo 
ies) 





| 4 71.8386 0.974768 
5 97.2413 0.865449 


6 


7 161.199 8.948593E-02 








135.027 


TRANSDUCER: 111/C30:HARD POINT/CESIUM UNIT Axis: +Z 
PEAK FREQUENCY (Hz) AMPLITUDE (g) 
1 7.81651 2.05755 | 
2 24.9373 9.57841 
3 60.0392 0.904517 


0.232790 





Table 5. Data Derived From Figure 7, 





IV. DESIGN TOOLS 
A. I-DEAS SOFTWARE 

I-DEAS [Ref. 4] is a product of the Structural Dynamics Research Corporation (SDRC) 
I-DEAS is an integrated package of mechanical engineering software tools. Among the 
applications available with the software are solid modeling design, drafting, finite element 
simulation, dynamics, post-processing and testing. 

In the vibrational analysis of the TOPAZ II power system, the finite element model 
solution, post-processing, and dynamics applications are the most useful. The faded solution 
application has the capability of creating a "solution set" which describes the type of solution 
(statics, dynamics, etc), the type of output to store, and other options. Using this application, the 
TOPAZ II finite element model was analyzed under both linear statics and dynamics. The 

| restraints used were a clamp located on each of the three feet, restraining the model against all six 
possible rigid body motions. The model was solved in order to obtain the natural frequencies of 
the structure. This solution also gave the displacements and stresses for each mode. These 
results were listed, and displayed two-dimensionally and color-coded, using the post-processing 
application. The modes were animated to determine each mode shape. Figure 8 is a snapshot of 
mode one. 

After the model was solved, the dynamics application was used to apply the vibration sine 
values of acceleration vs. frequency to the model. These loads were applied through the clamped 
feet. A viscous damping ratio of 0.02, determined from the previous modal test, was applied and 
the model response obtained and graphed. From these graphs the theoretical natural frequencies 


were obtained. [Ref. 4] 
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VALUE OPTION: ACTUAL 


1. OOE+03 


9. OLE+02 


8. OLE+02 


7. OLE+02 


6. OLE+02 


5. OLE+02 


4. OOE+02 


3. OGE+02 





2. QOE+02 


1. 00E+02 


0. OOF+00M © 


Figure 8. Mode 1, Displacement 
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V. FINITE ELEMENT MODELING AND ANALYSIS 
A. FINITE ELEMENT ANALYSIS - THEORY 

The finite element method provides a formulation which can use digital computer 
automation for the analysis of irregular systems. To accomplish this, the method sees a complex 
structure as an assembly of finite elements, where each element is part of a continuous structural 
member. By requiring that the displacements be computed and the internal forces in balance at 
points shared by several elements, called nodes, the entire structure reacts as a single entity. The 
finite element method expresses the displacement at any point of the continuous element in terms 
of a finite number of displacements at the nodal points multiplied by given interpolation functions. 
(Ref. 5] 

The advantage of the finite element method over other methods is that the equations of 
motion for the system can be derived by first deriving the equations of motion for a typical finite 
element and then assembling the individual elements' equations of motion. The motion at any 
point inside the element is obtained by means of interpolation, where the interpolation functions 
are the same for every element. Complicated structures can be modeled by simpler one or two 
dimensional elements, resulting in less computation time. [Ref. 6] 

The I-DEAS software utilizes the displacement based finite element method to solve finite 
element models. To derive the equilibrium equations, the principal of virtual work is used. The 
principle of virtual work is basically a statement of the static equilibrium of a mechanical system 
and is defined as follows: the work performed by the applied forces through infinitesimal virtual 


displacements compatible with the system constraints is zero. [Ref. 5] 
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The principle can be used to calculate the portion of static equilibrium of a system. For a 


conservative system, the equation is 


—— Caan N 
6W = Y F-8F = 23 — 6x + apt + ae = 0 (1) 


where 
8W = virtual work 

F, = the applied force 

ér, = virtual displacement vector of the 

dV = virtual potential 

8x; = virtual displacement in the X direction 
dy, = virtual displacement in the Y direction 
dzi = virtual; displacement in the Z direction 


As utilized by I-DEAS, the equation is as follows 


yt = (7st 4 FT of FT 
ee (o)dv = (U,) (F,) fe Ga aa Ga) (2) 


where 
( o ) = virtual displacement 
(© )= strain due to the virtual displacement 
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(€) = actual stress at current time t due to applied loads 

(F,) = nodal forces 

(f,) = surface forces 

(f, ) = body forces 

This equation states that the work done by point loads, surface loads, and body forces are 

equal to the internal strain energy due to the virtual displacement when integrated over the 
volume. Given a set of nodal values, the forces and displacements can be evaluated anywhere on 
the body. This is accomplished through the introduction of (N), continuous element interpolation 
functions, and (B), the oe of (N). Cancelling the virtual displacement vector from above, 


leaves 


= [BY (oa, =(F)+ ¥ fATapas, « 3 fayaar, 


elem V, elem S, elem V, (3) 


A global set of equations are formed from the integrations over each element. When considering 
i 


the linear case, this equation is reduced to 


[K](u) = #) + &,) + ®) (4) 


~ where 
(u ) = displacement 


{ K ] = linear stiffness matrix 
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The linear stiffness matrix is equal to 


[kK] = ¥ [By (De Iav, 


elem V7, (8) 
where 
|D] = matrix of elastic constants 
To determine strain and stress, the following equations are used 
(e) = [B]@) (6) 
(0) = [D](e) (7) 


A finite element model can be made of more than one type of element. The idea is to 
model the deflection of the structure, not its appearance. The FEM is the idealization of the entire 
structural problem. This includes the elements node locations, both the physical and material 


properties, boundary conditions and loads. [Ref. 4] 
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I. Modal Analysis 

The FEM was tested using I-DEAS software. The restraints used were clamps on each of 
the three feet. The model was solved, using the parameters from Appendix A, to determine its 
natural frequencies. Because the model used ideal clamps, these natural frequencies were 
expected to be higher than those found experimentally. No losses were modeled, also 
contributing to higher theoretical natural frequencies. These frequencies are listed in Table 6. 

When the FEM was solved, post processing revealed the results shown in Table 6. Using 
I-DEAS, modes up through 200 Hz were animated to determine each mode shape present. The 
upper bound of 200 Hz was chosen because that was as far as the experimental testing went. The 
mode shapes demonstrate how the structure moves when a particular natural frequency is excited. 
Figures 8 - 12 are snapshots of five of these mode shapes. Each snapshot gives the mode number, 
corresponding frequency, and color-coded scale relating to either stress or displacement. 
Displacement values are displayed in Figures 8, 9, 10, and 11, and stress values in Figure 12. 
Either parameter may be selected. 

The numerical stress and displacement values on these pictures are significant as relative 
quantities only. Using Figure 9 as an example, it is known that any "red" section will deflect ten 
times as far as any "royal blue" section of TOPAZ. Ifthe figure were displaying relative stress, 
the same would hold true for values of stress. 

When solving for natural frequencies, the -DEAS software is solving an eigenvalue 
problem. The shapes of the functions are known as the eigenfunctions. Because the equations 
being solved are homogeneous, the shape of the function can be deternmined uniquely but the 


amplitude of the functions are arbitrary. Each natural frequency has a corresponding 
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THEORETICAL DATA 


(Hz) AMPLITUDE | AMPLITUDE 


ANIMATED 
MODE SHAPE 


mode 2 beam, +/- Y 





mode 2 beam, +/- X 





7 [90s | szamem [6 | tatomei-¥, op 2 | 
: bottom 2% 

9 7 bottom only, radial +/- Z 
10 structure mode 3, XZ plane 


YZ deflection of reactor 
support 


81.94865 


99.2141 


4.3E+02 


3.55E+02 


torsion, frame 











150.7446 


22 





1.90E+02 
3.30E+02 
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13. | 115.4492 16 XY plane, deflection 2 
lower bays, opposite 

14 | 115.6353 mode 3, ZX plane 

15 120.2439 torsion in lower 2 bays 

ig | 127.1178 | 2.69B+02 

19 18 torsion about Z, top 
opposite bottom half 

20 130.4791 1.47E+02 


mode 3 ZX 
mode 3 ZX 





MODE 


FREQ 
(Hz) 





THEORETICAL DATA, CONT. 








DISPLACED | RELATIVE 
AMPLITUDE | AMPLITUDE 


26 177.8734 1.28E+02 
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34 


36 
3) 





178.1774 





ae eee 


182.0637 


183.776 
185.9247 


189.7822 


190.003 
193.3796 


35 193.6082 1.26E+02 ee i 


199.7701 1.18E+02 


1.15E+02 31-32 
1.57E+02 22-23 





8.96E+01 


cae 
Tasco [2 

; 28 
30 
3 


212.4125 1.10E+02 


Table 6. Theoretical Results. 
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ANIMATED 
MODE SHAPE 


torsion of bottom bay 
mode 3 ZY 
mode 4 ZX 
mode 4 ZY 
torsion, structure 
mode 3 ZX 
torsion, lower 2 bays 
torsion, structure 
torsion, structure 
axial, torsion, structure 


mode 3 ZY, torsion lower 
half 


axial, torsion 





VALUE OPTION: ACTUAL 


3. 38H+02 


3. 04E+02 
2. 70E+02 
2. 36B+02 
2. O38+02 
1. 69B+02 
1. 35B+02 
1. OLE+02 
6. 75E+01 
3. 38B+01 
0. OOE+00 


bo 


Figure 9. Mode 3, Displacement 
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VALUE OPTION: ACTUAL 


1. Q0E+03 


—_ . 00B+02 

= 
8. 00E+02 
7. OOB+02 
6. 00B+02 
5. ODE+02 
4. QOE+02 
3. OOB+02 
. OOE+02 
1. OOE+02 
Q. COE+00 


oO 


to 


Figure 10. Mode 6, Displacement 
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VALUE OPTION: ACTUAL 


1. O3B+02 


9. 29E+01 
ie 
8. 26E+01 
7. 22E+01 
6. 19E+01 
5. 16H+01 
4. 138+01 
-10E+01 
2. O6E+01 
1. O3E+01 
0. COE+00 


) 





bo 


Figure 11. Mode 38, Displacement 
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VALUE OPTION: ACTUAL 
SHELL SURFACE: TOP 


4 

Pf t 1 o * 
a eet ha 2. 238+06 
i a 1. 98E+06 
7 a 1. 13B+06 
ANA 1. 49B+06 
é 1. 24E+06 
9. 918+05 
T.448+05 
4. 96B+05 
2.498405 
1. 41B+03 


Figure 12. Mode 39, Stress 
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eigenfunction, or natural mode. Because the problem is homogeneous, any constant multiplier of 
a natural mode represents the same natural mode. Therefore, amplitudes are relative values only. 
[Ref. 6] 

Table 6 also lists the relative amplitudes of the displaced mode shapes. This is important 
to determine, when excited, which mode shapes will deflect the structure the greatest amount. 
This information also gives a feeling of where the maximum areas of stress will occur. A 
description of the animated mode shapes are listed in the last column of the Table. All axis 
references in this table are to the FEM coordinate system. Again, these may be used to predict 


the structural deflection when the corresponding frequencies are excited. 
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2. Forced Response 

In order to determine the accuracy of the FEM, it had to be tested in the same manner as 
was the actual structure. Table 7 gives the sine-vibration input that was applied to the actual 
TOPAZ unit. This same input was applied to the FEM. The sine vibration is shown graphically in 
Figure 6. This function was applied to the FEM through each of the clamped feet. A damping 


ratio of 0.02 was applied to all elements of the FEM, corresponding to the damping ratio of the 


actual structure. 














SINE VIBRATION 


Table 7. Applied Sine Vibration [Ref. 6]. 


Frequency ( Hz ) 




















The sine-vibration test was applied twice, once in the lateral direction and once in 
longitudinally, as was done in the experimental testing. The same lateral direction (Z of the actual 
unit) was used for testing of the FEM. 


When post processing the theoretical results, data was taken from locations on the FEM 
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which corresponded to accelerometer Jocations on the actual structure. As with the experimental 
data, the theoretical results were graphed. Each of these graphs were then compared to each of 
the experimental graphs obtained from the same location on TOPAZ. An example of one of the 
theoretical granhs may be seen in Figure 13. This graph came from node 92 of the FEM, which 
corresponds ¢o accelerometer number 10 of the actual TOPAZ structure. This accelerometer was 
mounted on the reactor top plenum. (See Table 3) As explained in Chapter IL, this comparison is 
necessary to determine the accuracy of the FEM. The results of this analysis are explained in 
Chapter V. 

In Figure 13, under the graph, the ID number of 3021 represents this particular function in 
I-DEAS. Each function was assigned a different number automatically. The functions were then 
plotted. The data from this figure were obtained from the X direction of the FEM. The Y axis 
labeled 'modulus' represents a non-dimensional number. The theoretical amplitude values were 
not necessary in determining which natural frequencies were excited. As with the experimental 
data, prominent peaks were tabulated from the graphs (See Figure 13 and the corresponding data 
listed in Table 8). The phase plots, such as the one in Figure 13, were useful in identifying which 
peaks were significant. When the phase changes 180 degrees, a natural frequency of the system 
has been crossed. As the frequency values increase, the phase plot becomes less angular and 
phase changes less distinct. This is because lower natural frequencies are also contributing to 
phase changes. 'X' marks are indicated on Figure 13 to show which peaks were tabulated. 

The remainder of the theoretical graphs and tabulated data may be found in Appendices B 


and C, respectively. 
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Figure 13. Sample Theoretical Plot. 
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PEAK FREQUENCY (Hz) 
1 10.0 
2 44.0 


i 6 150.0 
Sr 
7 176.0 











Table 8. Data Derived From Figure 13. 
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VI. ANALYSIS AND CONCLUSIONS 
A. ANALYSIS OF EXPERIMENTAL AND THEORETICAL RESULTS 

The theoretical and experimental graphs of Appendices B and D, respectively, were 
analyzed the same way. Prominent peaks, indicating natural frequencies of the system, were 
marked on the graphs. The raw data were then examined so that the frequencies at which these 
peaks occurred could be exactly identified. This procedure was followed for both the lateral and 
longitudinal tests of both the theoretical and experimental data. The lateral and longitudinal test 
results were analyzed separately. This was done because different natural frequencies would be 
excited during the two tests. 

At this point, there were four sets of tabulated results. These were the lateral and 
longitudinal results from the theoretical and experimental tests. Each of these sets were 
individually studied to determine which frequencies had peaks commonly occurring. The 
experimental lateral and longitudinal frequencies are listed in Table 9. The theoretical lateral and 
longitudinal results were combined and listed as the theoretically predicted natural frequencies of 
‘tte system. (See Table 9) 

As expected, the experimentally determined natural thea snise were lower than those 
theoretically predicted. As mentioned earlier, this is in part due to the ideal restraints used by the 
computer model. Also, no losses were simulated on the FEM, but losses would occur in an actual 
system. Some of the theoretically predicted frequencies were not found in the experimental data. 
This was due, in part, to the accelerometer placements on TOPAZ. They may not have been 
located properly to register each natural frequency. Some frequencies may have been damped out 


by losses in the system. This accounts for the "missing" experimental frequencies in Table 9. 
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NATURAL 


FREQ (Hz) 
1 9.2949 
2 9.4301 
3 27.6152 
4 40.57694 
44.111 
6 56.7768 
7 59.0425 
8 59.05082 


81.85867 
81.94865 










99.2141 
13 115.4492 








115.6353 

120.2439 
125.5667 
126.487 








130.483 
150.6242 
150.7446 
170.473 














THEORETICAL 


THEORETICAL DATA 


SINE-Z 
ACTUAL | ACTUAL 
(Hz) 

7.26 7.6 
8.8 
28 28 

z 
44 46 
53 $5 
58 56 
59 


101 
112 





121 
124 











ANIMATED MODE 
SHAPE 


beam, structure, +/- Y 
beam, structure, +/- X 
torsion, about Z 
mode 2 beam, +/- Y 
mode 2 beam, +/- X 
axial +/- Z 
bottom +/- Y, top +/- Y 


bottom +/- X 


9 63.30842 64 60 bottom only, radial +/- Z 
10 


structure mode 3, XZ plane 


YZ deflection of reactor 
support 





torsion, frame 


XY plane, deflection 2 
lower bays, opposite 


mode 3, ZX plane 
torsion in lower 2 bays 
axial. structure 


mode 3 beam, XY plane 


18 127.1178 mode 3 XZ 

19 129.8433 torsion about Z, top 
opposite bottom half 

20 130.479] 131 mode 3 ZY 


mode 3 ZX 
mode 3 ZY 
mode 3 ZX 
mode 3 ZX 






























THEORETICAL DATA, CONT. 


MODE | THEORETICAL | SINE-X 
NATURAL = | ACTUAL 
FREQ (Hz) (Hz) 


174.9176 171 torsion of bottom bay 


Panes [| is | moe zx | 
| 1820037 sf | torsion, structure 
Tesco [| _torn over2e— 
af ws7ea | | i ——ton te 
eo 
esse [Tin se 


193.6082 mode 3 ZY, torsion lower 
ee 


| 199.7701 7701 | axial, torsion | torsion 


re eae ae ae 


Table 9. Theoretical and Experimental Result Tabulation. 



















ANIMATED MODE 
SHAPES 
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Figure 14. Theoretical and Experimental Result Comparison. 
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The theoretically predicted natural frequencies were then plotted against the 
experimentally found frequencies. This graph is shown in Figure 14. For 70% of the predicted 
modes, there exists a corresponding experimentally determined natural frequency. Modes where 
experimental frequencies were not found were most likely a result of inadequate accelerometer 
placement. In a one to one comparison, theoretically predicted and experimentally found natural 


frequencies were within 4.6 Hz of each other, with an average difference of 1.8 Hz. 


B. CONCLUSIONS 

The finite element method is a very useful tool for structural analysis. Complex structures 
can be modeled, have forces applied, and be tested for failure in a number of ways, using only 
computer software. This provides a substantial cost savings when used to complement prototype 
or actual unit testing. The I-DEAS software proved invaluable for the analysis work presented in 
this thesis. 

In order to obtain reliable useful data from finite element modeling, the model must 
accurately represent how the actual structure will deflect under loading. One method of checking 
the accuracy of the model is to solve it for its natural frequencies. If the system modes predicted 
by the model resemble the actual system modes, the model is considered reliable. Through the 
testing completed for this thesis, the results show the two sets of natural frequencies, theoretical 
and experimental, to be very close, and therefore the model is an accurate representation of the 


TOPAZ structure. 
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C. RECOMMENDATIONS FOR FUTURE STUDY 

Further study of the finite element model test data should be done to locate the highest 
amplitude response, thereby identifying the frequencies which correspond to the worst case 
loading. As mentioned previously, individual components such as the cesium block or starter unit 
should be analyzed to see if they exceeded the limits of the test. If the power system must be 
launched in an upright position to conform to the shape of U.S. launch vehicles, structural 
analysis of the frame must be done to see if it will withstand launch loads. From previous testing 
in Russia, the interface between the shield and reactor can be a problem source. Further research 
needs to be completed in this area. 

When data were analyzed in this thesis, only relative values of stress and displacement 
were obtained. Further analysis needs to be done to predict the actual stress and displacement 


that will occur during a typical launch sequence. 
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APPENDIX A. 


This Appendix contains the parameters used in the I-DEAS [Ref. 4] software to solve the 


finite element model. 


Type of solution: 


BOUNDARY CONDITION SET 





Normal Mode Dynamics - SV1 
Restraint Set (Pinned Feet) 
Constraint Set (Base Constraints) 
Temperature Set 


OUTPUT SELECTION 


Mode Shapes 
Stress 


OPTIONS 


Method 

Mass Lumping 

Enforced Motion 

Number of Rigid Body Modes 


SOLUTION CONTROL 


Number of Flexible Modes 
Number of Iteration Vectors 
Maximum Iterations 
Frequency Accuracy (digits) 
Start Vector Creation 
Frequency Shift Method 
Frequency Shift in Hertz 


ET 


Normal Mode Dynamics 


Selected 
Selected 
Selected 
Room Temperature 


Store/List 
Store/List 


Solution No Restart 
Selected 

Selected 

0 


43 

58 

10 

3 

Automatic 

Shift and Number of Modes 
0 


A minimum of 10 iterations were required to converge on a solution. 
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APPENDIX B. 


This Appendix contains 54 theoretical graphs obtained using the I-DEAS [Ref. 4] 
software. The FEM was solved after applying a sine vibration input. Results were taken from 
nodes of interest, corresponding to accelerometer locations on the actual unit. These results were 
then plotted and the graphs are contained herein. 

Pages 54-80 contain graphs from the longitudinal-sine input. Pages 81-107 contain graphs 


from the lateral-sine input. 
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APPENDIX C. 


This Appendix contains the theoretical results obtained from analysis of 25 of the graphs 
from Appendix B. Only the graphs corresponding to experimental testing accelerometer locations 
were analyzed. Prominent peaks were located and the corresponding frequencies tabulated. A 
table was made corresponding to each of the 25 graphs. 

Pages 110-113 contain results from the longitudinal-sine input. Pages 114-117 contain 


results from the lateral-sine input. 


109 


De aa a 
X-AXIS SINE (LONGITUDINAL) THEORETICAL RESULTS 


NODE 92-Y 
FREQUENCY (Hz) 
178.0 





NODE 92-X 
PEAK FREQUENCY (Hz) AMPLITUDE 
1 10.0 





44.0 
58.0 
82.0 
116.0 
150.0 
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|} 


176.0 


NODE 97 - X 
FREQUENCY (Hz) AMPLITUDE 
1 10.0 


; 


58.0 
4 82.0 


6 150.0 


ve 176.0 
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NODE 97-Z 


FREQUENCY (Hz) AMPLITUDE 
10.0 
44.0 
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PEAK FREQUENCY (Hz) AMPLITUDE 
ees ee (ee 


ise ee 
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NODE 55 - X 
FREQUENCY (Hz) AMPLITUDE 
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Z-AXIS SINE (LATERAL) THEORETICAL RESULTS 
NODE 75 - Z 
PEAK FREQUENCY (Hz) AMPLITUDE 
ee 
2 44.0 - 
3 58.0 ; 
4 176.0 : 








NODE 79 -Z 


176.0 





NODE 97 - X 


FREQUENCY (Hz) 










AMPLITUDE 














114 






ee 
pee ee 











NODE 5 - X 
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a ae eee ea 
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NODE 55 -Z 


PEAK FREQUENCY (Hz) AMPLITUDE 














1 8.0 - 
2 42.0 2 
3 56.0 - 
4 64.0 e 
5 100.0 “ 


NODE 29 - X 
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APPENDIX D. 


This Appendix contains experimental graphs obtained from testing of the actual TOPAZ 
unit. Data were taken from accelerometers placed at various locations on the unit, manipualted 
by a MATLAB routine, and then plotted. 

Pages 120-143 contain graphs from the longitudinal-sine input. Pages 144-167 contain 


graphs from the lateral-sine input. 
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Amplitude, g 





Transducer: A1:RESP.;INBD.SPIDER FX.LEG Axis: +X 
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Transducer: A2:CTL.;TBL.FX.NEAR INB.SPDR. Axis: +X 









Amplitude, g 
3, 














10° 
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Amplitude, g 





Transducer: 108/C01:REAC. LEG BRACKET/+Y Axis: -X(20) 














2 
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10° 10 
Frequency ~ Hz 
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Amplitude, g 








Transducer: 110/C02:REAC.TOP PLENUM/+Y Axis: +Y 











Frequency ~ Hz 
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Amplitude, g 


Transducer: 110/C03:REAC.TOP PLENUM/+Y Axis: -X 








Frequency ~ Hz 
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Transducer: 107/C05:REAC. LEG BRACKET/-Z Axis: -X(20) 











Amplitude, g 
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Frequency ~ Hz 





Amplitude, g 





Transducer: 109/CO9:REACT. TOP PLENUM/-Z Axis: -X 
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Amplitude, g 


Transducer: 109/C12:REAC.TOP PLENUM/-Z Axis: -Z 











Frequency ~ Hz 
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Amplitude, g 


Transducer: 101/C13:-Z LEG @ BASE Axis: -X 











Frequency ~ Hz 
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Amplitude, g 


Transducer: |01/C14:-Z LEG @ BASE Axis: -Z 











Frequency ~ Hz 
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Amplitude, g 





Transducer: 101/C15:-Z LEG @ BASE Axis: -Y 
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10 10 10 
Frequency ~ Hz 
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Amplitude, g 
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Transducer: 102/C16:LEG/+Y AND +Z/BASE Axis: -X 


























10 10 10° 
Frequency ~ Hz 
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Amplitude, g 





Transducer: 102/C17:LEG/+Y AND +Z/BASE Axis: +Y 











i?) 2 3 
10 10 10 10 
Frequency ~ Hz 
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Amplitude, g 


Transducer: 106/C19:LEG/+Y AND +Z/RAD. TOP Axis: -X 
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10° 10 
, Frequency ~ Hz 
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Amplitude, g 





Transducer: |06/C20:LEG/+Y AND +Z/RAD. TOP Axis: +¥(30) 
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10° 10 10 
Frequency ~ Hz 
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Transducer: 104/C24:BOTTOM COLL.RAD./+Y Axis: +X 


Amplitude, g 
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Amplitude, g 


Transducer: 103/C27:BOTTOM COLL.RAD./-Z Axis: +X 

















Frequency ~ Hz 
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Amplitude, g 





Transducer: |111/C29: HARD POINT/CESIUM UNIT Axis: -Y 
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Amplitude, g 


Transducer: 111/C30:HARD POINT/CESIUM UNIT Axis: -Z 

















Frequency ~ Hz 
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Amplitude, g 


Transducer: 105/C31:JOINT ON LEG Axis: -X 














Frequency ~ Hz 
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Transducer: 112/C34:STARTUP UNIT FRAME Axis: -X 
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Amplitude, g 


Transducer: 112/C35:STARTUP UNIT FRAME Axis: -Z(60) 
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Amplitude, g 


Transducer: 112/C36:STARTUP UNIT FRAME Axis: -Y(60) 

















Frequency ~ Hz 
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Amplitude, g 








Transducer: A1:SINE CTL.1;INBOARD LEG(N.) Axis: +Z 








10 10° 
Frequency ~ Hz 
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Transducer: A2:SINE CTL.2; OUTBOARD LEG(E.) Axis: +Z 
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Amplitude, g 


Transducer: 108/C04:REACT. LEG BRAKET/+Y Axis: +Z 

















Frequency ~ Hz 
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Transducer: 107/CO8:REACT. LEG BRAK./-Z Axis: -Z(15) 
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Amplitude, g 
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Transducer: 109/C12:REACT. TOP PLENUM/-Z Axis: -Z 





Frequency ~ Hz 
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Amplitude, g 


Transducer: 101/C13:-Z LEG @ BASE Axis: -X 











Frequency ~ Hz 
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Amplitude, g 





Transducer: 101/C14:-Z LEG @ BASE Axis: -Z 
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Amplitude, g 


Transducer: 101/C15:-Z LEG @ BASE Axis: -Y 














Frequency ~ Hz 
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Amplitude, g 


Transducer: 102/C17:LEG/+Y AND +Z/BASE Axis: +Y 














Frequency ~ Hz 
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Amplitude, g 





Transducer: 102/C18:LEG/+Y AND +Z/BASE Axis: -Z 
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Transducer: |06/C19:LEG/+Y AND +Z/RAD. TOP Axis: -X 
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Amplitude, g 
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106/C20:LEG/+Y AND +Z/RAD. TOP Axis: +Z(30) 
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Amplitude, g 








Transducer: 104/C23:BOTT. COLL. RAD./+Y Axis: +Y 








Frequency ~ Hz 
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Amplitude, g 


Transducer: 111/C28:HARD POINT/CESIUM UNIT Axis: -X 
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Amplitude, g 





Transducer: !111/C29:HARD POINT/CESIUM UNIT Axis: -Y 









Frequency ~ Hz 
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Transducer: 105/C31:JOINT ON LEG Axis: -X 
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Transducer: 105/C33:JOINT ON LEG Axis: -Z 
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Amplitude, g 


Transducer: 112/C34:STARTUP UNIT FRAME Axis: -X 
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Amplitude, g 


Transducer: 112/C35:STARTUP UNIT FRAME Axis: -Z(60) 
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APPENDIX E. 


This Appendix contains the experimental data obtained from analysis of the graphs of 
Appendix D. Prominent peaks were located and the corresponding frequencies tabulated. A table 
was made corresponding to each graph. 

Pages 170-175 contain results from the longitudinal sine-vibration test. Pages 176-183 


contain results from the lateral sine-vibration test. 
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TOPAZ 4 X-AXIS SINE (LONGITUDINAL) EXPERIMENTAL DATA 
TRANSDUCER: A1:RESP.;INBD.SPIDER FX.LEG Axis:+X 


MODE FREQUENCY (Hz) AMPLITUDE (g) 














TRANSDUCER: A2:CTL.;TBL.FX.NEAR INB.SPDR. Axis: +X 


MODE FREQUENCY (Hz) AMPLITUDE (g) 


TRANSDUCER: 108/C01:REAC. LEG BRACKET/+Y Axis: -X(20) 


TRANSDUCER: 110/C03:REAC.TOP PLENUM/+Y Axis: -X 


MODE FREQUENCY (Hz) AMPLITUDE (g) 
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TRANSDUCER: 107/C0S: REAC. LEG BRACKET/-Z Axis: -X(20) 


TRANSDUCER: 109/C09: REACT. TOP PLENUM/-Z Axis: -X - 


MODE FREQUENCY (Hz) AMPLITUDE (g) 


TRANSDUCER: 109/C12: REAC.TOP PLENUM/-Z Axis: -Z 



























0.933774 


TRANSDUCER: 101/C14:-Z LEG @ BASE Axis: -Z 


2.36357 







1.00910 
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TRANSDUCER: 101/C15:-Z LEG @ BASE Axis: -Y 


FREQUENCY Gi) | __ANPLITUDE 
l 


7.26411 4.428413E-02 





15.1196 0.217778 


39.0735 1.24001 


102/C16: LEG/+Y AND +Z/BASE Axis: -X 


7.26411 0.956400 


2 38.0587 1.68669 
103.69 0.884210 


TRANSDUCER: 102/C17:LEG/+Y AND +Z/BASE Axis: +Y 


ort cr 
4 L 


















TRANSDUCER: 
AMPLITUDE (g) 



































154.734 0.250545 


TRANSDUCER: 106/C19:LEG/+Y AND +Z/RAD. TOP Axis: -X 
FREQUENCY (Hz) AMPLITUDE (g) 








. 7.26411 0.947660 


39,5012 5.00209 











107.438 0.670597 
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TRANSDUCER: 104/C24: BOTTOM COLL.RAD./+Y Axis: +X 






















MODE FREQUENCY (Hz) ' AMPLITUDE (g) 
[— 1 7.26411 0.916813 
2 38.4494 6.19911 















1.61763 





112.567 








TRANSDUCER: 103/C27: BOTTOM COLL.RAD./-Z Axis: +X 
FREQUENCY (Hz) AMPLITUDE (g) 

7.26411 0.957385 

4.70091 



















38.4494 















109.031 





1.38565 









FREQUENCY 


TRANSDUCER: 111/C29: HARD POINT/CESIUM UNIT Axis: -Y 
AMPLITUDE (g) 































1.32583 








40.9289 2.26666 















2.83179 


1.94570 
0.913474 


58.1368 
84.0117 
















123.977 
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Wie ee eens, = Peet See A Eh i aioaea eg pee ot 
TRANSDUCER: 111/C30: HARD POINT/CESIUM UNIT Axis: -Z 


FREQUENCY (Hz) AMPLITUDE (g) 


8.67318 3.010930E-02 


69.7589 0.597402 


113.349 0.901246 

























TRANSDUCER: 105/C31: JOINT ON LEG Axis: -X 


















TRANSDUCER: 112/C34:STARTUP UNIT FRAME Axis: -X 


ee 
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TRANSDUCER: 112/C35:STARTUP UNIT FRAME Axis: -Z(60) 


MODE FREQUENCY (Hz) AMPLITUDE (g) 


25.2811 0.312611 





33,8624 0.434271 
39.7120 0.862081 
52.9984 0.437670 
71.5108 0.342384 
101.493 0.246964 
150.347 0.162090 
181.722 0.234352 


TRANSDUCER: 112/C36:STARTUP UNIT FRAME Axis: -¥(60) 
MODE FREQUENCY (Hz) AMPLITUDE (g) 
7.25590 0.114642 
34.0583 3.08350 
38.9664 6.18548 


101.493 0.546350 





172.147 0.555203 








Z-AXIS SINE (LATERAL) EXPERIMENTAL DATA 


TRANSDUCER: A1:SINE CTL.1;INBOARD LEG(N.) Axis: +Z 
FREQUENCY (i) | __ AMPLITUDE @ 
1 7.5904 0.750334 


2 28.4307 0.895611 


3 169.748 0.844949 


TRANSDUCER: 













A2:SINE CTL.2; OUTBOARD LEG(E.) Axis: +Z 


FREQUENCY (Hz) AMPLITUDE (g) 
7.5904 0.777016 


AMPLITUDE (@) 
ets I 


TRANSDUCER: 107/C08: REACT. LEG BRAK./-Z Axis: -Z(15) 


AMPLITUDE @) 
| 
| 


89.8642 0.268446 


ae, ey 163.492 5.619752E-02 


































~ 
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TRANSDUCER: 109/C09: REACT. TOP PLENUM/-Z Axis: -X 


3 
4 98.3362 





er 
8.513646E-03 


TRANSDUCER: 109/C11:REACT. TOP PLENUM/-Z, Axis: + 
FREQUENCY (Hz) AMPLITUDE (g) 
8.76736 1.95501 








MODE 
15.3975 9.63723 1E-02 
15.7615 8.896597E-02 
28.3918 2.16219 


54.1565 0.603319 


59.8494 0.212064 
69.5546 0.141099 


84.2027 0.188112 


111.574 2.786758E-02 
10 ~ 118.594 2.285344E-02. 
11 143.565 2.197792E-02 








12 190.176 1.028353E-02 
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TRANSDUCER: 109/C12: REACT. TOP PLENUM/-Z Axis: -Z 


MODE FREQUENCY (Hz) 


29.4789 
53.8231 
100.092 
133.096 
175.661 





ttt 





AMPLITUDE (g) 


9.40651 4.61358 


3.97788 
1.36627 
0.138697 
5.688959E-02 
1.350727E-02 





TRANSDUCER: 101/C13:-Z LEG @ BASE Axis: -X 


MODE FREQUENCY (Hz) 

7.59040 
24.5878 
28.0407 
32.1489 
46.3899 
56.3866 
68.8677 
81.9021 
85.1832 
120.925 


—_ 











1] 138.978 
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TRANSDUCER: 101/C14:-Z LEG @ BASE Axis: -Z 


AMPLITUDE (g) 
0.548643 
9.492084E-02 
0.114665 
8, 86981 8E-02 
4.471105E-02 
5.088181E-02 
4.105160E-02 
1.659305E-02 
2.000362E-02 
0.140088 
0.265800 
0.110414 


















TRANSDUCER: 101/C15:-Z LEG @ BASE Axis: -Y 


PLTODE 

































AMPLITUDE (g) 
0.502636 
0.625009 
0.728032 
0.683575 





: 





26.8651 
131.239 1.16256 


TRANSDUCER: 106/C19: LEG/+Y AND +Z/RAD. TOP Axis: -X 
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TRANSDUCER: —_106/C20: LEG/ +Y AND +Z/RAD. TOP Axis: +Z(30) 
MODE FREQUENCY (Hz) AMPLITUDE (g) 
7.75082 1.56896 
24.1758 2.22697 
28.8881 2.56046 
55.9942 0.349077 


96.9556 0.466277 





136.480 0.899289 


TRANSDUCER: 104/C23:BOTT. COLL. RAD./ +Y Axis: +Y 
FREQUENCY (Hz) AMPLITUDE (g) 
9.86712 0.105801 





MODE 


27.943 2.34252 
48.6237 4.5289 
69.3226 2.20957 
88.3314 3.07370 
130.777 2.95019 





154.658 2.88126 


















TRANSDUCER: 103/C25:BOTT. COLL. RAD./-Z Axis: -Z 


FREQUENCY di 


MODE 
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TRANSDUCER: 111/C28: HARD POINT/CESIUM UNIT Axis: -X 


MODE AMPLIRIDE 


TRANSDUCER: 111/C29: HARD POINT/CESIUM UNIT Axis: -Y 


MODE FREQUENCY (Hz) AMPLITUDE (g) 
8.83778 1.23716 
38 


















































TRANSDUCER: 111/C30: HARD POINT/CESIUM UNIT Axis: +Z 


MODE FREQUENCY (Hz) AMPLITUDE (g) 
2 24.9373 9.57841 
3 
5 





























a ern 
a 
Patina 
SO 
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TRANSDUCER: 105/C31: JOINT ON LEG Axis: -X 


MODE FREQUENCY (Hz) AMPLITUDE (g) 
7.76132 0.410090 


51.9497 0.122066 
73.0626 0.256050 















TRANSDUCER: BSL eesOINe ON LEG Axis: -Z 


PLITUDE 
a ST 


TRANSDUCER: 112/C34:STARTUP UNIT FRAME Axis: -X 


AMPLITUDE 


TRANSDUCER: 112/C35:STARTUP UNIT FRAME Axis: -Z(60) 


a 
ee 
ee 
oe er | ee 
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0.831606 



















TRANSDUCER: 














112/C36:STARTUP UNIT FRAME 





Axis: -Y(60) 





AMPLITUDE (g) 



































1 7.97612 1.13312 
2 24.2569 1.43728 
28.3918 1.85856 
59.4898 0.371952 








108.367 





139.493 
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0.362733 





0.220601 
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